Large enhancement of the thermopower in Na^Co02 at high Na doping 
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Research on the oxide perovskites has uncov- 
ered electronic properties that are strikingly en- 
hanced compared with those in conventional met- 
als. Examples are the high critical tempera- 
tures of the cuprate superconductors and the 
colossal magnetoresistance in the manganites. 
The conducting layered cobaltate NaxCo02 dis- 
plays several interesting electronic phases as x 
is varied [J, 0, including water-induced su- 
perconductivity [4\ and an insulating state [3| 
that is destroyed by field [5j. Initial measure- 
ments [3| showed that, in the as-grown compo- 
sition, NaxCo02 displays moderately large ther- 
mopower S and conductivity a. However, the 
prospects for thermoelectric cooling applications 
faded when the figure of merit Z was found to be 
small at this composition (0.6< x <0.7). Here we 
report that, in the poorly-explored high-doping 
region x >0.75, S undergoes an even steeper en- 
hancement. At the critical doping Xp ^ 0.85, Z 
(at 80 K) reaches values ^40 times larger than in 
the as-grown crystals. We discuss prospects for 
low-temperature thermoelectric applications. 

In the large- a; region of NaxCo02 {x >0.75), progress 
has been hampered by difficulties in growing single crys- 
tals as well as by phase-separation effects which appear 
above Xp. Powder neutron diffraction |fi] has revealed 
subtle shifts in the Na ions as x is increased above 0.75. 
The doping interval 0.75 < x < Xp is a homogeneous 
phase H2, in which the thickness t of the C0O2 layers 
undergoes a steep increase 0. Above Xp, the neutron re- 
sults indicate phase separation (the mixed-phase region is 
labelled H2+H3). Muon spin rotation nuclear mag- 
netic resonance (NMR) and susceptibility 0, 
experiments also suggest mixed phases in highly-doped 
samples. Finally, the limiting phase H3 (x ^ 1) is in- 
ferred from NMR and magnetic susceptibility to be a 
non-magnetic insulator with a slight admixture of a resid- 
ual metallic phase 0, Q| . 

In a thermopower experiment, the current density 
J' = (tE produced by the £^-field is cancelled by the 
thermoelectric current J = a(— VT) driven by the ap- 
plied gradient — VT (a the thermoelectric or Peltier con- 
ductivity). Whereas the thermopower S — a/ a is the 
quantity usually reported, we have found that a pro- 
vides a more incisive probe for sorting out the results in 
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the region H2 + H3. Results from 10 crystals (of nominal 
size 400x700 /im) with x ranging from 0.71 to ~1 reveal 
that a rises steeply to a peak at the threshold doping 
Xp, with a nearly constant (as shown in Fig. ^ samples 
are labeled 1-10 in order of increasing x). In the mixed 
region above Xp both quantities decrease rapidly. 

Crystals of NaxCo02 grow as either bilayer or tri-layer 
structures (2 and 3 C0O2 layers per unit cell, repec- 
tively) 0. The in-plane transport quantities are very 
similar in the 2 structures. However, at large x, the 3- 
layer crystals (Samples 2*, 9*, 10*) tend to lock to the 
commensurate dopings x = j and ^1, which provide 
valuable checks. 

Figure n displays curves of the observed p vs. T (tem- 
perature) in the 3-layer (Panel a) and 2-layer (b) cobal- 
tate. Beginning with the lowest curves (Samples 1-4 in 
Panel b), we note that the proffies of p{T) are metallic 
with the residual resistivity ratio RRR ~ 20, comparable 
to the as-grown crystals (0.6 < a; < 0.7). Significantly, 
in Samples 2*-6 (regions H2 and slightly beyond Xp), 
the room-temperature values p(300) lie in the range 1-2 
mi7cm. As x crosses Xp into H2 + H3 (Samples 5-8), 
p(300) rises steeply to ^--^240 mfJcm (in 10*), which far 
exceeds the Mott limit. Surprisingly, the curves retain 
the "metallic" profile even as a; ^ 1~. This implies that, 
for X > Xp, the carriers strongly segregate into a fraction 
/ of the layers. The observed p{T) reflects the intrinsic 
resistivity of these conducting layers inflated by the geo- 
metric factor 1//. Moreover, the persistence of metallic 
behavior to 4 K implies that the layers remain electrically 
connected down to low T. Thermal activation of carri- 
ers over low barriers, frequently seen in phase-separated 
systems, is not observed. 

Curves of the thermopower are displayed in Fig. 
The curve with the smallest S (Sample 1) is typical of 
as-grown crystals P, Q . Across the H2 region (Samples 
2*, 3, 4), the S-T curves show an increasingly negative 
curvature that grows into a broad peak of magnitude 240 
pV/K at ^^130 K in 4. In the mixed-phase region, the 
peak value rises further to 300-350 pY/K (Samples 5-8), 
before settling to 200 pV/K as a; ^ 1 (Samples 9*, 10*). 
These large S values far exceed any reasonable extrapola- 
tion of the Sommerfeld expression S — {ks / e){kBT / ep) 
(ks is Boltzmann's constant and e the electron charge). 
To match the observed dS/dT ~ 3 pY /K^ at low T, we 
would need the Fermi energy to equal ~30 K, which 
is unphysically small. The unusually large S, coexist- 
ing with low resistivities (^ 100 /if2cm in Samples 4-6), 
challenges our understanding of how strong correlation 
enhances the Peltier effect. Approaches include the use 
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FIG. 1: In-plane resistivity p vs. T in the 3-layer (Panel 
a) and 2-layer (Panel b) cobaltate in log-linear scale. The 
samples, numbered in order of increasing a:, are: Sample 1 
{x ~ 0.71), 2* (0.75), 3 (0.80), 4 (0.85), 5 (0.88), 6 (0.89), 
7 (0.96), 8 (0.97), 9* (0.99) and 10* (1.0). Asterisks indi- 
cate the 3-layer crystals (Panel a). From x — 0.71 to Xp ~ 
0.85, p does not change significantly (Samples 1-5). Inside 
the mixed-phase region H2 + H3, however, p increases rapidly 
(6-10*). The persistent metallic profile inside H2 + H3 sug- 
gests that conducting layers are embedded in an insulating 
matrix. Scaling of the curves of p are not as satisfactory as 
for a (see below). Contacts with contact-resistance 2-10 
were attached by lightly abrading the crystal surface and at- 
taching Au wires with Ag paint (Du Pont 4922N). The size 
of contact pads (see Supplement) leads to a total uncertainty 
in the absolute values of p of ±15%. 



of the Heikes formula 0, 111113 • 

In Fig. we display plots of the figure of merit 

Z = S^/pK for some of the samples (k is reported in Sup- 
plementary Information). Near the threshold Xp (Sam- 
ples 4, 5), Z rises to a prominent maximum below 100 K. 
As mentioned, the peak value of Z at 80 K is ~40 times 
larger than that in the as-grown composition (compare 5 
with 1). We discuss Z below. The dimensionless figure 
of merit ZT is plotted in Supplementary Information. 

The Peltier conductivity a shows nominally similar 
profiles (Fig. In Samples 1 and 2*, a displays broad 
peaks at 100 and 80 K, respectively, but remains fairly 
small. As we enter the H2 region, however, a increases 
rapidly, displaying a sharply peaked profile in Sample 4. 
In simple metals, a ^ {k'^/e){T/eF)cr decreases linearly 
to zero as T ^ 0. At high T, a tends to saturate to 
values in the range 1-10 A/mK. Hence, the pronounced 
peak and large values of a reported here {amax 135 



400 



300 



50 100 150 200 250 300 



100 




2.0 



1.5 



1.0 



0.5 



0.0 







a 








J A 
















■■■■■■■■i 




5 


; CBT 


b 


. / 4 






- f/3. 








1 •♦♦♦♦♦♦♦^ 





50 100 150 200 250 300 
T(K) 



FIG. 2; The in-plane thermopower S (Panel a) and the figure 
of merit Z (Panel b) in NaxCo02. In Sarnple 1 (Panel a), S 
is very similar to that of Terasaki et al. As 2; increases 
into phase H2 (2*-4), the profile of S develops an increasing 
bulge near 130 K that grows rapidly to peak values of 200-250 
/iV/K. In the mixed-phase region, S further increases to 300- 
350 AtV/K (Samples 5-8), before settling down to 228 p.V /K 
in the limit a; — > 1 (Samples 9*). A striking pattern is that 
the S-T profiles (all hole-type) are nominally similar in shape 
in the mixed-phase region. The exception is Sample 10* in 
which / ~ 1/200 (see Supplementary Information). Panel 
b shows curves oi Z = / pn, with k measured separately 
(not shown). As x increases from 0.71 (Sample 1), the peak 
value of Z increases steeply to 1.8xl0~^ in Sample 5. 

At higher x, the peak value of Z falls rapidly because of the 
sharp increase in p. At 80 K, Z in Sample 5 is ~40 times 
larger than that in Sample 1. The dashed line labeled CBT 
is Z reported [13 for CsBi4Te6. 



A/mK in Sample 4) are strikingly anomalous, and likely 
a consequence of strong correlation. 

By tracking T^ax (the peak temperature of a), we see 
that a vs. T in 1-4 {x < Xp) evolves differently from a 
in 5-10* (in H2 + H3). In the first g roup, Tjyidx decreases 
from 110 K to 50 K (in 1 and 4, respectively). The contin- 
uous change in Tmax precludes a simple scaling analysis. 
By contrast, in the second group, Tmax remains fixed at 
~ 50 K. Moreover, by multiplying each curve by a scale 
factor p, we may match it to the curve of 4 (Fig. ^p). 
The exception is Sample 6 which has a cusp at 36 K (of 
unknown origin) that ruins the scaling. Hence, as we 
traverse the homogeneous region 772, the peak value of 
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FIG. 3: Curves of the Peltier conductivity a vs. T (Panel 

a) and scaling of the a-T curves in selected samples (Panel 

b) . As a; traverses the region H2 (Samples 1-4) the peak in 
a(T) systematically shifts from 100 K to 50 K. Because of the 
shift these curves cannot be scaled together. However, in the 
region H-i, + Hz (5-10), the peak in a remains at 50 K. Panel 
b shows that the curves of 5-10 match the template in 4 when 
multiplied by the scale factor p {p = 1.3, 3.9, 5.4, 38 and 200 
in Samples 5, 7, 8, 9*, and 10*, respectively). Sample 6 has a 
cusp of unknown origin that makes it an exception. Sample 
10* is discussed in Supplement. 
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FIG. 4: Variation of a (50 K) and a (300 K) vs. x in 
NaxCo02 at large x (Sample labels shown). As x traverses 
H2, Q:(50) rises sharply to a maximum at Xp in correlation 
with the increase in the C0O2 layer-thickness t 0. Because 
cr(300) is nearly constant, the peak in Q!(50) leads to a peak 
in Z as well. In H2 + H3, a and a (plotted vs. Xm calculated 
from p using Eq. decrease linearly as Xm ^ 1~, implying 
phase-separation of holes into conducting layers. The phase 

9*, 10*) X is 
1. The X values 

in Samples 1 and 4 were measured by x-ray diffraction (lines 
are guides to the eye). Error bars in Xm reflect the uncertainty 
in measuring p and S (see Supplementary Information). 



Hz is insulating. In the trilayer crystals (2 
locked to the commensurate values | and 



a increases rapidly from 20 A/mK in 2* to 135 A/mK 
in 4, but the form of a vs. T also changes continuously. 
However, once we cross the threshold Xp into the mixed 
phase, the profile is locked to that of 4. 

The simplest explanation of the scaling is that, in 
H2 + H^, the carriers segregate into continuous conduct- 
ing layers embedded in a nearly insulating matrix. The 
conducting layers correspond to doping Xp whereas the 
insulating matrix is at x = 1.0. The mean doping Xm is 
then given by 

Xm= fxp + {l- f), {Xp<X^<l) (1) 

where / is the fraction of layers with doping Xp. Addi- 
tivity of the Peltier currents implies that the observed a 
is given by 

a(T) = fap{T) + (1 - /)ai(T), (2) 

where subscripts p and 1 refer to quantities evaluated at 
X = Xp and 1, respectively. 

The scaling behavior in Fig. implies that the second 
term in Eq. |21 is negligible if / is not too small. The 
Peltier conductivity ai(r), which has a maximum value 
of ~ 0.85 A/mK, is observable only when / becomes very 
small, as in the case of Sample 10* (see Supplementary 
Information). Away from this limit, we may take / ~ 1/p 
throughout the mixed-phase region; this allows us to find 



Xjn using Eq. ^ For e.g., in Sample 9*, the scaling of a 
shows that 1 in 40 of the layers is conducting {xm = 0.996 
instead of 1.00), which is consistent with the observed 
p(30G) in Fig. In addition, we have employed other 
checks on the calibration. The trilayer crystals, with x = 
I and ~1, provide valuable calibration points. Moreover, 
X in Samples 1 and 4 was determined from the c-axis 
lattice parameter measured by x-ray diffraction. These 
checks lend support for our calibration. 

With the inferred values of x and Xm, we next plot the 
variation of a at 50 K (circles) and a (triangles) at 300 
K in the phase diagram at large doping (Fig. QJ- From 
X — 0.71 to Xp, a(50) rises by a factor of 14. The modest 
changes in (t(300) (triangles) and k (see Supplementary 
Information) result in a substantial increase in Z. We 
emphasize that these increases occur within the homo- 
geneous phase H2 and are possibly correlated with the 
increase in layer thickness t Clearly, they are not a 
consequence of the phase separation that onsets at Xp. 
Both quantities decrease linearly above Xp. 

The 40-fold enhancement of Z realized at 80 K im- 
proves greatly the prospects for thermoelectric applica- 
tions. The peak value in Fig. 03 (Z - 1.8 x 10"^ 
K^^) is among the highest known for a /loZe-type mate- 
rial below 100 K. The promising material CsBi4Te6 also 
displays a large Z at 100 K In comparison, Z in 
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NaxCo02 peaks at low temperatures whereas the curve 
for CsBi4Te6 (curve labeled CBT) falls steeply. The al- 
loy Bii_,ySby has long been known to display an 
even larger Z at low T (at 80 K, Z = 7 x IQ-^ K^i 
for y ^ 0.12). However, it is electron-like. To realize 
its advantages in a Peltier device operating below 100 
K, we need to pair Bi-Sb with a hole-type material with 
comparabe Z, but the 35-year search has turned up no 
viable candidates. We calculate that an optimized de- 
vice with NaxCo02 (x — Xp) paired with Bii_ySby has 
a device figure-of-merit Znp — 2.5 x 10^^ K^^ at 80 K 
[Znp = [Sn - SpY{^pnK.n + ^ PpUp)'"^ where subscripts 
n and p index the 2 materials]. This value is higher 
than any reported to date. Further enhancement of Z 
in NaxCo02 seems possible if we degrade k. 
Methods 

Sodium metal (0.5 g) and tetrahydrofuran (THF) (30 



mL) in a pyrex tube were warmed at 100 °C in an oil 
bath 0. Benzophenone (2 g) was added and the solu- 
tion heated until it turned blue from the formation of 
the benzophenone ketyl radical anion. As-grown crystals 
of Nao.75Co02 were placed along a pipette-like tube and 
immersed in the blue solution of the sodium radical an- 
ions. The pyrex tube was quickly capped and heated for 
4 days at 100 °C. Finally, the crystals were extracted and 
washed with THF and ethanol to eliminate trace sodium. 
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Supplementary Information 

Breakdown of scaling in Sample 10* 

We provide more details on the measurements and anal- 
ysis of the thermopower S and resistivity p in the phase- 
separated region H2 + H3, especially in the extreme case 
X ^ 1^ represented by Sample 10*. 
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FIG. 1: Comparison of the scaled curves of pa{T) in Sample 
10* with those in the batch Samples 4-9* (the scale factor p 
= 200 in Sample 10*). Above 100 K, the scaled curve matches 
those in 4-9*. Below 100 K, the additional contribution which 
we identify with ai in the x = 1 phase becomes prominent. 
The peak value of ai is 0.85 A/mK. This additional term 
causes the thermopower S{T) in 10* to be distinct in profile 
from Samples 4-9* (see Fig. 2a of the text). 



with Sp = ap/ffp. 

In Samples 5-9*, in which the curves of a vs. T may 
be scaled to that in Sample 4, the term (1 — f)c(i must 
be negligible compared with fax^ in Eq. O (Sample 6 
is excluded). Further, if the second term in Eq. |5l is 
neghgible (i.e. / >0.05), all the profiles of S{T) should 
be identical to Sp (Sample 4). We find that, apart from a 
20% spread in magnitude, the curves of S{T) are indeed 
closely similar in Samples 4-9*. 
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As noted in the text, the simplest interpretation of the 
scaling behavior in the curves of the Peltier conductivity 
a vs. T is to assume that the hole carriers strongly seg- 
regate into parallel layers throughout each crystal. The 
layers must remain robustly connected electrically even 
at 4 K (otherwise we would not have the metallic profile 
in p with the large residual resistivity ratios RRR^20). 
Additivity of the currents implies that the observed a 
and a are given by 

a(T) - /ap(T) + (l-/)ai(T), (3) 
<j{T) = }ap{T) + {\- f)a^{T), (4) 

where subscripts p and 1 refer to quantities evaluated at 
X = Xp and 1, respectively. 

The similarity of the p-T profiles in Fig. 1 of the text 
implies that all of J is carried by the conducting layers. 
Setting fji — 0, we have for the observed thermopower 

5(r) = 5,(T) + ii— (5) 



FIG. 2: Curves of the in-plane thermal conductivity k vs. 
T in Samples 3-8 of NaxCo02 (Panel a) and in Sample 10* 
(Panel b). The magnitude of k = 8-10 W/mK above 50 K 
in Samples 3-7 is characteristic of doped layered perovskites. 
In Sample 10* (a: ~ 1), k attains much larger values, possibly 
refiecting the pristine nature of the Na ordering. It is of the 
same order as crystals in the charge-ordered insulating state 
{x = 0.5). 

This spread reflects difficulties in establishing the ab- 
solute values of S and p accurately. The size of the con- 
tact pads (50-100 pm) relative to that of the crystals 
(400-700 pva on a side) leads to uncertainties of ±10% 
in estimating the voltage lead separation. This is exacer- 
bated by the large electrical anisotropy pd p where pc is 
the c-axis resistivity, which causes a small admixture of 
Pc into the measured "in-plane" resistance despite delib- 
erate care at optimal contact pad placement. We suspect 
that the latter factor is the main factor that spoils ex- 
act scaling of the p-T profiles. Altogether, we estimate a 
combined uncertainty of ±15% in the absolute values of 



6 



0.15 




FIG. 3: Plots of the dimensionless figure-of-merit ZT in Sam- 
ples 1-8 (except 2). 

p. Within this uncertainty, we see that, in Samples 4-9*, 
the profiles of S vs. T remain strikingly similar despite 
the 50-fold increase in p in these samples. 

The prominent exception is Sample 10*. As shown 
in Fig. n the scaled curve pa{T) matches those in the 
other samples above 100 K. The value of the scale factor 
p ~200 inferred from this plot is also consistent with 
the resistivity data {p{3Q0) in Sample 10* ~ 240 mficm 
nominally equals p(300) in Sample 4 multiplied by p). 
Below 100 K, however, a further contribution, which we 
identify with ai(T) in Eq. O becomes dominant. We 
note that the actual peak value of ai is ^ 0.85 A/mK, 
which is unresolvable in Fig. 3a of the text. While this 
is a very small value, it can make the second term in Eq. 



O comparable to the first if / becomes very small. This 
is the case in Sample 10* near 50 K, for which fap ~ 

0. 68 A/mK. Hence, in this extreme limit, we estimate 
that 1 in 200 of the layers carries the applied current to 
produce a "metallic" resistivity profile resembling that in 
Sample 4 (apart from the scale factor). In principle, the 
observed thermopower and a should also match that of 
4. However, in Sample 10*, the term in ai dominates 
that of the conducting layers below 100 K and ruins the 
scaling behavior. 

Thermal conductivity 
For completeness, we show the in-plane thermal con- 
ductivity measured in some of the samples in Fig. S2. 
These curves were used in computing the figure of merit 
Z = S'^/pK, with K the in-plane thermal conductivity. 
Throughout the region k has the modest value 10 
W/mK, quite typical of layered perovskites (Fig. S2a). 
However, in the limit x = 1, k becomes very large (Fig. 
S2b), attaining values in the charge-ordered state^ at x 
= 0.5. 

Dimensionless Figure of Merit ZT 

Some authors prefer to use the dimensionless figure-of- 
merit ZT in place of Z. For comparison, we have replot- 
ted the curves of Z vs. T in Fig. 2b of the text as ZT 
vs T in Fig. S3. An interesting feature of these curves 
is that ZT attains a prominent maximum at tempera- 
tures between 60 and 100 K in the layered cobaltate. By 
comparison, in high-performance thermoelectric materi- 
als based on doped Bi2Te3 the maximum in ZT occurs in 
the interval 200-300 K (see for e.g. Chung et al. The 
rapid fall of ZT in the Bi2Te3 family reflects the semi- 
metallic nature of the bands. In the cobaltates, however, 
the peaking of ZT at much lower T is a consequence 
of the very narrow bands and the magnetic correlations 
that are present because of strong interaction between 
carriers. 
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